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We investigate statistics of large-scale structures from large-eddy simulation (LES) of 
turbulent channel flow at friction Reynolds numbers Rct = 2k and 200k. To properly 
capture the behaviour of large-scale structures, the channel length is chosen to be 96 
times the channel half-height. In agreement with experiments, these large-scale structures 
are found to give rise to an apparent amplitude modulation of the underlying small-scale 
fluctuations. This effect is explained in terms of the phase relationship between the large- 
and small-scale activity. The shape of the dominant large-scale structure is investigated 
by conditional averages based on the large-scale velocity, determined using a filter width 
equal to the channel half-height. The conditioned field demonstrates coherence on a scale 
of several times the filter width, and the small-scale-large-scale relative phase difference 
increases away from the wall, passing through 7t/2 in the overlap region of the mean 
velocity before approaching tt further from the wall. We also found that, near the wall, 
the convection velocity of the large-scales departs slightly, but unequivocally, from the 
mean velocity. 



1. Introduction 

Re cent studies ( Kim fc Adrianlll999 [ Morrison. McKeon. Jiang fc Smitsll2004l:l iGuala. Hommema fc AdrianI 



•ian 
de^ 



2006 : lMontv. Stewart. Wilhams fc Chongl2007l:lH utchins fc Ma.rusic'2007&"a':'Mathis. Hutchins fc Marusi^ 
2009) have confirmed earlier observations ( Favre. G aviglio & Dumas 1967; Kovaszna v. Kibcns fc Blackwelde: 
1970 ) of very long large-scale structures in the wall region of boundary layers, channels 
and pipes. These structures are marked by streamwi se-elongated, alternating low- and 
high- moment um, meandering st reaks, with length IPS { Kim fc Adrianll999t Morrison et al. 



2004) ■ 8-16 J (iGuala et a/ll20G6D 256 (iMontv et ^l2007h. 6^ (iHutchins fc Marusi(j20076h " 
20^ (iHutchins fc Marusidl2007a[ ) and width 0.3-0.5^ (jMathis et aZ.ll2009l). where J is the 
boundary layer thickness, channel half- height or pipe diameter. See Monty et al. ( 2009() 



for a description of the differences between the characteristics of these large structures in 
the different canonical flows. The bursting period of these structures, 6S/U, where U{z) 
is the mean velocity, was already noted some decades ago, along with the long tails of 
the streamwise velocity auto-correlations, see review by Cant we J ( 198l[ ). The dynamical 



significance of these large-scale structures can be seen in a scale decomposition of relative 
energy content, as measured by the premultiplied one-dimensional longitudinal spectrum 
KxEuu plotted against the log streamwise wavelength, logA^^, where A^; ~ 2n/ (equal 
area under the curve implies equal energy contribution). For a boundary layer at fric- 
tion Reynolds number i?e^ = 7.3k ( Hutchins fc Marusic 2007&() . the signature of these 
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structures are related to the outer peak in KxEuu found at {z/6, \x/5) = (0.06, 6). It has 
been proposed that the wall-norr nal location of this peak is located at the middle of the 
log layer, z+ oc Rey^ or Rel^^ (jMathis et all |2009| ) (the choice of seal ing depends on 
whether the lower limit of the log law is Reynolds number dependent), where z is the 
height from the wall and the superscript + indicates scaling in wall units: the friction 
velocity it,, and kinematic viscosity v. However the scaling remains somewhat ambiguous. 



The large-scale structures were found (jBandvopadhvav fc Hussainlll984t iMathis et al 



2009( ) to modulate the amplitudes of superimposed small-scale fluctuations. To test this 



idea, these authors first split the streamwise velocity into large- and small-scale com- 
ponents via a filter at Xx/S = U{z)/{fS), and then used either a filtered and rectified 
small-scale signal or the Hilbert-transform to determine the envelope for the small-scale 
fiuctuations, finally forming the correlation coefficient between the large-scale fluctua- 
tions and the low-pass filtered envelope of the small-scale fiuctuations. They found that, 
near the wall, large-scale high-speed regions carry intense superimposed small-scale fluc- 
tuations, but this correlation is reversed above a height that decreases in outer units with 
Rcr- We shall attempt to reproduce these features presently. 

The footprint of structures centred far from the wall provides an obvious challenge 
in terms of determining appropriate convection velocities across the range of turbulent 
scales, with particular importance for obtaining the correct wavenumber spectra from 
temporal frequency spectra obtained by, for example, hot-wire anemometry. It has been 
known for some t ime that convection v elocities deviate from the local mean in the near- 
wall region, (e.g. Krogstad et al. 1998h . The common practice is to use Taylor's frozen- 
turbulence hypothesis to map from the frequency to the wavenumber domain, that is to 
use the assumption that all structures at a given wall distance z convect at the same 
scale-independent mean velocity U(z). It w as shown from a particle image velocimetry 
(PIV) experiment ( Dennis &: Nickel^ 2008[) that this is indeed good approximation at 
z/S = 0.16 for a Reg = 4.7 k boundary layer, at least for scales smaller than their field of 
view, 3.2(5 in space and 6.3S/U in time. However note that this wall- normal distance is 
sufficiently far from the wall that it is beyond the large-scale energy peak, such that any 
convection velocity questions are likely insignificant because of the low shear in the outer 
region. With a field of view larger than 20(5 x 20S/U and height down to z/5 — 0.049, we 
revisit the question of whether the footprint of the large-scale structures, having centres 
further from the wall, still convect at the local mean velocity near the wall. 



25(5 (Montv et al. 


2007) 


(Chune & PullinI 


2009li. 



LES since its cost depends only on the number of 'large eddies', which, for a channel, is 
Reynolds number independent. In contrast, the fully resolved direct numerical simulation 
(DNS) is prohibitively expensive. For re ference, the most ambit ious DNS of a channel 
flow to date is the Rer = 2k simulation ( Hovas Sz Jimene"3 2006 ) in an Lx/S = Stt « 25 
channel, where Lx is the streamwise length; a DNS investigation at higher Rer and 
larger Lx of these large-scale structures is not yet possible. Of course, the use of LES 
comes at the cost of subgrid-scale (SGS) modelling, wall modelling and numerical errors, 
but LES is much faster (hours for simulation, minutes for post-processing) than DNS 
and experiments; we hope that a controlled application of the present LES combined 
with experience in the subject may shed some light on the physics of these large-scale 
structures. 

Details of the simulations are given in §[2] and discussion of observations arc found in 
§|3] before we conclude in §|4l 



LES of large-scale structures in long channel flow 



3 



Run Rer L^/S Ly/5 h+ A^/S Atu^/S Ny Nt TUc/L^ 

Gl 2 k 95 7.9 15 0.17 0.006 576 48 48 72000 110 
H3 200 k 96 8.0 750 0.083 0.002 1152 96 96 15800 12 
Gib 2 k 95 7.9 15 0.17 0.006 576 48 48 22000 34 

Table 1. LES parameters for long channel flows. Channel-transit times based on 
data-recording period T and centreline velocity Uc- 
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Figure 1. Comparison between DNS data of iHovas fc Jimenej (j2006l l and 
present Re = 2 k channel flow LES, run Gib (table [!}. Spectra at 
z/S = 0.049, 0.090, 0.17, 0.26 <t^> z+ = 98, 180, 350, 510 (ordered in decreasing energy). 



2. Simulation details 

As full details of t he LE S, including the numerical method and SGS model, are given 
by IChung fc Pullinl ( 20091 ). we only highlight the important points here. We solve the 
filtered Navier-Sto kes equations for the LES velocity field u using the stretched-spiral 
vortex SGS model (|Misra fc PullinlllQQTtlVoelkl et alhoO(i\ . To circumvent the inhibitive 
cost of resolving the n ear- wall region (IChapmanlll979l) . z < ho, we use a wall model 
( Chung fc Puiiinll2009l) that supplies ofF-wall slip- velocity boundary conditions at ho to 
the interior LES, operating in ho < z < 2S — ho, where z = 0,25 locates the walls. 
Presently, we fix ho = 0.18 A^, and the slip velocity is calculated using a wall model 
comprising 1) an evolution equation for the wall shear stress derived from assuming 
local inner scaling for the streamwisc momentum equation and 2) an extended form of 
the stretched-vortex SGS model that provides a local log relation along with a dynamic 
estimate for the local Karman constant. 

The parameters for the three LES runs are given in table [T] The grid is uniform. 
Ax = Ay = 4:Az, throughout the simulation domain. To capture the physics of long 
large-scale structures, we use long a channel, L^/S w 96, and the statistics are taken over 
T channel-transi t times. 

It was shown (jChung fc Pullinll2009l ) that the LES-prcdicted statistics from the Re = 
2 k case, includi ng means, turbulent int ensities and spectra, compare reasonably well 
with the DNS of iHovas fc Jimenez ( 2006 ). We show the root-mean-square (r.m.s.) of the 
streamwise velocity fluctuations and the LES-resolved spectra in figure [TJ The points 



4 



D. Chung and B. J. McKeon 



in[lja) correspond to actual discretization points. Even though the total (subgrid plus 
resolved) r.ni.s. is within 90% of the DNS result, its spectra plotted in energy-content 
form overpredicts the DNS spectra by about 20%. (Note, however, that the comparison 
between experimental channel flow data and DNS has revealed a similar trend, with the 
experimental premultiplied spectra at large waveleng ths in the overlap laye r being up 
to order 10% larger than the equivalent DNS values ( Montv &: Chong |2009| ).) As such, 
the results presented here should be viewed as approximate, despite capturing energy 
at significantly larger wavelengths. On the other hand, the physics reported here can be 
seen as robust features of wall turbulence if they are also observed elsewhere. We note 
that the peak value, (kxEuu)^ « 1 at z/S = 0.090, X/S ~ 6 (figure [T]) is within the 
range of the peak va l ues from boundary layer spectra for Rer = 1.0-7.3 k (see figure 9 of 
Hutchins fc Marusic ( 2007qI )). We interpret the LES results as a model of the real flow. 



and emphasize that Rct = 200k is far out of the reach of current DNS approaches. 

In order to compute correlations based on temporal averaging, a numerical rake in runs 
Gl and H3, fixed in streamwise-spanwise location, is set up to record the LES velocity 
u and its modelled subgrid fluctuations T^x (= uu — uu) at the wall-normal locations 
z = UzAz {uz = 0,1,..., Nz) and time steps t = utAt {rit = 0,1, Nt — 1). Analogous 
correlations based on spatial averages are also computed from these runs from a snapshot 
in time. 

The only difference between runs Gl and Gib is the recorded data. For the latter, 
the three-dimensional data set, u{nxAx,y,nzAz,ntAt), is recorded at fixed y, for = 
0, 1, ... , 575, = 0, 1, . . . , 48 and = 0, 1, . . . , 21999, where A^/S = 0.17, Az/S = 0.041 
and AtUr/S = 0.006. 

Presently, x, y and z respectively denote the streamwise, spanwise and wall- normal 
directions; the velocity components, u, v and w, are defined accordingly. 



3. Discussion of observations 

We present observations of the LES velocity fields, with emphasis on the large scales. 

3.1. Convection velocities from spatio-temporal spectra 

We begin by using the somewhat unique combination of spatial and temporal data avail- 
able in this study to investigate the validity of Taylor's hypothesis. Given the autocorre- 
lation of the streamwise velocity fluctuations, 

R{p,t) = {u{x,t)u'{x + p,t + T)), (3.1) 

where p is the streamwise separation; r is the time delay; and u = U + u' such that 
{u') = 0, we define the spatio-temporal spectrum ^{k,uj) as the Fourier transform of 
R{p, t). That is, together they form the Fourier transform pair, given by 

^{k,u;) = j^J^^J ^R{p,r)e-'^-P-^'>dpdr, 

/OO pOO 
/ «'(K,a;)e'('^''-"^MKdw. 
-OO J —OO 

The wavenumber spectrum 0(/s:) and the frequency spectrum ^{oj) are both related to 
^ via 

/oo POO 
*(K,w)dw, $(w) = / «'(K,w)dft, (3.2) 
-CO 'J — oo 
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whence the mean-square of u fluctuations can be recovered from 

/ ^{K,Lu)dKduj = / $(w)dcj= / e{K)dK. 

-oo J — oo J —oo J ~oc 

These are even, $(a;) = $(— oj) and 6(k) = 6(— k). Since can be measured directly 
using hot wires, the one-dimensional longitudinal spectrum, 

e(K) EE (3.3) 

strictly defined in terms of a dispersion relation and associated group velocity 

ui'{k) = dw/dK, is often approximated using Taylor's frozen-turbulence hypothesis. Under 
this assumption, the dispersion relation can be formally written as ojt{k) = Uk and, from 
(|3.3p . Q{k) = U^{Uk). This is another interpretation to Taylor's hypothesis based on 
the Fourier decomposition: the group velocity ui'rp^K) = U and phase velocity lut/k = U 
of all eddies contributing to the turbulent kinetic energy at a particular wall-normal 
location are constants independent of wavenumber and equal to the mean velocity at 
that location, an approximation accurate for sufhciently small eddies. To obtain G from 
this strict interpretation of Taylor's hypothesis can be relaxed to account for energetic 
eddies which travel at U ±AU provided ^'(k, uj) is symmetric with respect to the uj ^ Uk 
line, that is 

^{K,ij) ^^{uj/U,Uk), (3.4) 
because p.4p then relates the two definitions in 



e(K) EE / ^{K,Uj)du}= ^{u}/U,Un)duj =U ^{K',UK)dK' = U^{Uk). 

J —oo J —oo J —oo 

For the purpose of obtaining wavenumber spectrum from frequency spectrum, a test of 
the validity of Taylor's hypothesis can be recast as a question of the symmetry of ^'(k, lS) 
with respect to the line lo = Uk. When this symmetry is broken, a different dispersion 
relation u}c{k) is necessary to relate 4> to 9. Since we have access to both $ and 9, we 
compute the loc{k) directly by finding the inverse to the monotonic function Kt 

poo poo 

uJc{k) ^ K^\Kx{k)), Kx{k)= 9(K')dK', Kt{uj)= <i>{oj')d^\ 



where Kx{0) = Kt{0) = (u'^)/2. The analogy to the convection velocity of individual 
eddies is complex. To first order, 0Jc{k) describes the apparent passing frequency of the 
most energetic eddies with wavenumber k, although this is an integral effect over the 
range of energetic spanwise wavenumbers. 

When computing 4* from the present LES simulation, a normalised Hann window in 
time, \/2/3 [1 — cos{2nnt / Mt)] {MtAt is the temporal window size), is applied to u before 
taking the discrete Fourier transform because the velocity is not periodic in time. From 
run Gib (table [ij where Nt = 22000 and Mt = 576, the spectrum is averaged across 
[22000/(576/2) — Ij =75 half-overlapping windows. No windowing is necessary in the 
periodic streamwise direction. 

We plot in figure [2] contours of the premultiplied spectrum, k^ki'^/u^ versus logA^; 
and \og\t, where A^: = 2n/Kx, \t = ^nj Kt and Kt = ujjU . When Taylor's hypothesis 
is valid, contours of n^Kt'^'/u^ should be symmetrical about the At = A^: line. Observe 
from figure [2] that Taylor's hypothesis is indeed a good approximation, except near the 
waU, z/5 = 0.041 and 0.083, and for the large scales, Xt/6, X^/S > 10. Note that the LES 
formulation docs not permit effective examination of smaller scales or locations closer to 
the wall. The dispersion relation line computed from At = XxUi/uic, appearing below 



6 D. Chung and B. J. McKeon 




Figure 2. Premultiplied spatio-temporal spectra of streamwise velocity fluctuations, 
KxK.t'^/u^ = 0.1,0.2,0.3,0.4, at various heights from Rcr = 2k channel flow LES, run Gib 
(table [1]): , dispersion relation computed from ^, At = X^uj/ujc, , Taylor's hypothe- 
sis \t = \x. 
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(to the right) of Taylor's hypothesis, At = XxUj/lot ~ A^: • 1, imphes that the phase 
velocity Uc = uJc/k is larger than the mean velocity U. 

Despite overwhelming similarities, the p hase yelocity Uc is not strictly the same as the 
convection velocity defined by (2.7) from IWillsl ( 1964I ). where it is related to the ridge 



of Recall from (|3.3p that by construction Uc{k) is defined such that the wavenum- 
ber spectrum can be recovered from the frequency spectrum $, that is = 
^{KUc{K,))d{KUc{K)) /dn. Put another way, it is the mapping from wavenumber space to 
frequency space such that the energy observed in wavenumber space at the wavenumber 
K, Q{k) dn, is equal to that observ ed in frequency spa.ce ^( KUc(K))d(KUc(K)) at frequency 
hUc{k) (see also discussion in Montv fc Chonel ( 20091 )). In practice, see figure [2l this 



UJ 



integr al-based definition often traces out the ridge of that is the definition in IWillsl 
( 19641 ). Because it can be measured readily, the ridge-based definition is often taken as 



the surrogat e for Uc- Energy transmission occurs at neither Uc nor the velocity given by 



Willsl (|1964[ ). but the group velocity d{KUc{n)) / dn, that is the velocity of the energy of 
wavepackets with wavenumber near k. Summarising, to obtain the wavenumber spectrum 
from the frequency spectrum, we require both the phase velocity and group velocity, see 
p.3p : the former appearing in the argument of $ as kUc{h) to rescale the frequency; and 
the latter appearing as a factor of $ as d{KUc{K)) / dn to rescale the rate of change in 
frequency. 

Following iDennis fc Nickeld (l2008l) . we can also test the validity of Taylor's hypothesis 



in physical space by plotting the autocorrelation R defined by (JSTTj). Taylor's hypothesis 
is valid, or more precisely there is a straightforward conversion from the temporal to 
the spatial domain, where the contours of R are symmetrical about the px = pt line, 
where pt = tU . Like figure[2l R in figure[3]shows an unequivocal departure from Taylor's 
hypothesis near the wall, z/5 = 0.041 and 0.083 and for the large scales, pt/S, px/5 > 10. 
There also appears to be a slight discrepancy of the opposite sign for the large scales 
at z/S ~ 0.5, which is more marked in the autocorrelation t han our presentation of the 
spectrum. The boundary layer PIV experiment performed bv lDennis fc Nickelsl (l2008[ ) at 



Ree w 4.7k reported that Taylor's hypothesis is still valid at the height z/6 = 0.16 for 
the field of view p/S < 0.29m/0.09m = 3.2 and tU/S < Is x 0.57ms-70.09m = 6.3. 
This is consistent with the present LES data since at z/5 = 0.17, figure [3l[c), R is indeed 
symmetrical about the px = Pt line, even up to very large scales pt, Px = 20(5. 

Consideration of spatio-temporal spectra permits some speculation about the convec- 
tion velocities of the energetic structures and the error associated with identifying the 
footprint of eddies with a particular streamwisc scale on the near-wall region from tem- 
poral data. The LES velocity fields unequivocally indicate that the most energetic large 
structures convect faster than the local mean velocity close to the wall, with the deviation 
growing close to the wall. Conversely, these large structures convect slower than the local 
mean velocity near the channel centre, figure [3l^e, /). This suggests that these eddies are 
"local" to a region in the overlap layer, in the sense that the mean velocity matches their 
convective velocity somewhere in the log region. We speculate that the location of this 
velocity matching corresponds to the location of the large-scale streamwise energy peak, 
which is consistent with the approximate magnitude of the difference between lo / lot and 
w/wc for large wavelength. This suggests that the departure from Taylor's hypothesis at 
the large scales should strengthen with increasing Reynolds number due to the increasing 
shear near the wall. This is an area of current experimental study. 

3.2. Large-scale-small-scale interaction 
With the differences between the spatial and the temporal decompositions in mind, we 
now describe the correlation that characterizes the interaction between the large-scales 
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Figure 3. Spatio-temporal correlations of streamwise velocity fluctuations, 
R/u^ — 0.3,0.6,0.9,1.2, at various heights from Rcr = 2k channel flow LES, run Gib 

(tableO: , dispersion relation computed from vf, pt = pxijj/ujc\ , Taylor's hypothesis 

pt = px- 
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and the small-scales. We first extract the large-scale fluctuations ul by applying a sliding- 
window top-hat time average, centred at t, to u: 

I ft+T/2 

UL{t) = - / u{t')<lt', (3.5) 

T Jt-T/2 

where r is the width of the sliding window. A low-pass filter, p.Sp . dampens fluctuations 
with frequencies higher than l/r. In spectral space, p.5|) is equivalent to a multiplication 
by the filter sin(Ta;/2)/(aj/2), where ui is the angular frequency. For clarity, we have 
suppressed the a;-dcpcndcncc of u in this part of the discussion since x is held constant. 

The small-scale fiuctuations are defined to be the remaining part of the motion, us = 
u — ul. The small-scale intensity can be measured by its local r.m.s., 

us{t)= i-J^ 4(0di'j . (3.6) 

Physically, us measures the local envelope or intensity of small-scale fluctuations. For 
example, if us is normally distributed, 95% of its amplitude is estimated to lie within 
2us. If us is not normal, us still measures the spread or envelope of us- In any case, it| 
appears in the equati ons governing ul, obtained b y applying the filter p.5|) to the Navier- 
Stokes equations (see Reynolds fc Hussain ( 19721 ) for a related two-scale decomposition). 



wh ich is another way to interpret us - Not e that equivalent approache s have been used 
by iBandvopadhvav fc Hussainl (|l984h and iGuala. Metzger fc McKeonI (|2009l) to obtain 



similar results in a range of flows including a laboratory turbulent boundary layer and 
the near-wall region of the near neutrally stable atmospheric surface layer, respectively. 

An elegant alter native to obtain the envelope of us is via the Hilbert transform 
( Mathis et ai 20091 ). and this approach has led to a signiflcant advance in understanding 



of the large-small scale interactions. However it is harder to relate the results to the 
governing equations of turbulence. 

When calculating an r.m.s. defined locally, (|3.6p . one has to contend with the inevitabil- 
ity that large-scales (low frequencies) have been aliased into the small-scale signature. 
Perhaps a better alternative is to use a tapered window in p.6p . alleviating some but 
not all of the aliasing. We have tried this and found some minor changes, but the gen- 
eral picture is unaltered, and so we decided to keep the simple definition in p.6p . The 
Hilbert transform bypasses this aliasing difficulty at the enveloping stage, but the issue 
reappears when one filters the envelope signal. We note that even if a perfect decompo- 
sition can be found, nature herself does not permit it, that is the two peaks in KxEuu 



(jHutchins fc Marusidl20076[ ) are never completely isolated, at least in Fourier space. 



In terms of LES quantities, we can write (|3.5|) and (|3.6p as 



1 



t + T/2 



UL{t) = - / u{t')<lt', (3.7a) 

T Jt-T/2 

us{t) - - / [ul{t') + Txx{t')\ At' , (3.76) 



T 



t-T/2 



where u is the resolved velocity; us = u—ul; and T^x is the modelled subgrid fluctuations 
associated with time scales smaller than the numerical discretization At- Using p.7|) . 
we now construct the normalised large-scale-small-scale correlation based on temporal 
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filtering: 



R (.') = (K - U){us - {us))) „^ 
{{uL-U)^)y-^{{us-{us)Y)^/^' ^ ' 



where the global or ensemble average is formally given by 



1 



T->oo T 



T/2 



{4>) = lim - / 0(t')dt'. (3.9) 



-T/2 



Note that the visual envelope of ws, e.g. 2us; does not affect Rr because the constant 
factor cancels out in the normalised correlation p.Sp . In other words, does not contain 
amplitude information; it does, however, contain phase information, since it is the cosine 
of the angle (or phase) between ul — U and us — (us), using the inner product ( ). In 
practice, we obtain Rt{z) by replacing the integrals with sums and ensuring the recording 
period T is much larger than the largest physical time scale in the flow, see table [T] 

The spatial counterpart to p.8|) . Rp{z), is defined analogously, with x and p respec- 
tively replacing t and r, while holding other variables constant. For Rp, the global average 
(|3.9p is replaced by an average over the wall-parallel plane with area L^Ly. The spatial 
correlations are calculated from runs Gl and 113 at one snapshot in time. 

The physical meaning of the correlations R = Rr, Rp are as follows. If large-scale 
higher-speed regions carry higher small-scale intensity (positively correlated, in phase), 
then i? « 1. Similarly, if large-scale higher-speed regions carry lower small-scale intensity 
(negatively correlated, re out of phase), then R « — 1. R sa can occur either if there is 
no correlation between the large and small scales, or if they are n/2 out of phase, which 
is physically the more likely option given the strong cor r elatio n for small and large z/5, 
as already demonstrated bv iBandvopadhvav fc HussainI ( 1984 ): Mathis et al. (l2009l ). 



Although the r.m.s. -based correlation coefficient p.Sp is d ifferent from i ts Hil bert- 
transform-based counterpart in the boundary layer study of Mathis et al. ( 2009[ ). we 



expect similar qualitative features if the large-scale-small-scale phase relationship is a 
universal aspect of wall-bounded flows, namely channels and boundary layers. As we are 
interested in large scales with sizes of the order 105, we expect that the LES will capture 
this statistic satisfactorily since the premise of an LES is to directly simulate the large 
scales; presently, there are 6 grid points per 5 in both wall-parallel directions and 24 grid 
points per 5 in the wall-normal direction. 

Figure [4] compares the correlations based on temporal filtering, Rr , and correlations 
based on spatial filtering, Rp, for Rcr = 2 k and Rcr ~ 200k and different values of r 
and p. Observe that near the wall, and us are positively correlated, up to i? « 0.4 
(the maximum in the domain we resolve, although note that the maximum value likely 
increases closer to the wall), but above a certain crossing height, z/S w 0.2 for Rcr — 2k 
and z/S ^ 0.11 for Rcr ~ 200 k, they are negatively correlated, down to i? ~ —0.4. The 
trend of decr e asing crossing height with increasing Reynolds number is also reported by 
Mathis et al. ( 20091 ) for the turbulent boundary layer, with z/S « 0.07 for Rcr = 2.8k 



and z/S 0.03 for Rer = 19 k. Presently, the correlations are largely independent of filter 
sizes 1.7 < tU/S,p/S < 12.3, although the deviation between Rp and Rr with the larger 
filter sizes from the one with p/S = 1.7 close to the wall is exacerbated in the spatial plots, 
as would be expected from the arguments concerning Taylor's hypothesis at the large 
scales in the preceding section. The kick-up of Rp relative to Rr for curves corresponding 
to p/S — 7.0, 12.3 in the vicinity of z/S = 0.45, see figurelUa, b) is presumably related 
to the convection ve locity effect d e monst rated in figure [Sj^e). A small sensitivity to filter 
sizes is reported by iMathis et all (|2009f l. Perhaps a precise quantitative comparison is 



impossible owing to the different envelope-extraction techniques and the different type 
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of wall-bounded flows. Note that a robust feature is that R increases slightly at the 
centre of the channel, but remains negative (figure U]). This increase is also reported in 
Mathis et al. ( 2009[ ) (and this sensitivity appears to be slightly enhanced in the case 



of the spatial correlation), but their increase is from negative correlations to positive 
correlations, a feature possibly related to the intermittent boundary layer thickness not 

present in channel flows. 

Although the large- scale-small- scale interaction was recently ( Mathis et aZ.ll2009h framed 



in terms of amplitude modulation, the results could also be discuss ed in terms of the rela- 

tive ph ase between the large and small scales, as originally posed bv lBandvopadhvaY &: Hussain 



(|1984[ ) and implied by the formulation of (|3.8I 



3.3. Conditionally averaged large-scale velocities and small-scale intensities 

To gain some insight into the structure of the large-scale coherent regions and the phase 
relationship between ul and us shown in figure IH we now turn our attention to condi- 
tionally averaged ul and us fields computed with the streamwise filter window p = 5. As 
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Figure 5. Conditionally averaged streamwise velocity fields in the (a, 6) spanwise-wal- 
1-normal and (c, d) streamwise-wall-normal planes, with (uJ[j4i)(Aa;) and {u'^\Ai){Ax), 
where Ai is the low-speed event ul{xx) — U < with probability P{Ai} = 0.53 at 
Xx/5 = (0,0,0.090): (a, c) filled contours of \{u+\Ai} - C/+| = 0.1,0.3,0.5,0.7, line contours of 
|(Sj|yli>-(5:J>| = 0.01; and (6, d) filled contours of Kwjl^i) - ("s)! = 0.01,0.02,0.03,0.04, line 
contours of \{u'^\Ai) — U^ \ = 0.1. Darker to lighter shades and dashed lines for negative values. 



Vectors represent in-plane velocity components. Data from snapshot of Rct 
LES, run Gib (table [l]). 



2 k channel flow 



seen m figureHand in lMathis etal\ (|2009t ). this phase relationship is relatively unaffected 
by the choice of p. 

Figure [5] shows ensemble averages, {ul\Ai){Ax) and {us\Ai){Ax), conditioned on the 
occurrence of a large-scale low-speed event at z/J = 0.090, Ai = ul{xx) — U < 0, where 
Aa; = X — Xx, Xy^ = (0,0,180) ^ Xx/S = (0,0,0.090), computed from a snapshot 
of LES run Gib (tabl e [H . The spanwise-wall- nor mal view, figure [5][a), is previously 
shown in figure 7(a) oflHutchins fc Mari^ (|20076h with the choice = 150 from the 
DNS data ofldd Alamo et "dl (|2004 at Rbt ~ Ik. For reference, the vectors in figure 
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[5fa, b) correspond to LES discretization points. The agreement between DNS and LES 
results is striking. The essential features of the DNS averages are also seen in the present 
LES averages, namely 1) a splatted low-speed region with minimum « —0.8 and 

width A.y/S « 0.4 centred on Xx flanked on both sides by weaker high-speed regions with 
maximum w 0.1; and 2) the in-plane large-scale swirl at A^/S w ±0.2, z/6 ~ 0.2. 

We report that P{Ai} = 0.53, implying that nearly identical figures, but with signs 
reversed, are seen when we condition on the large-scale high-speed event A'l ~ uj^^ > 
(the complement of ) because 

0^{<j,)^{^\A,)P{A,} + {cj,\A[)P{A[} ^ {^\A^)^~{cj,\A[), 

(equality holds if P{Ai} is exactly 1/2). Since the reversed picture, conditioned on A[ 
exists, neighbouring high-speed regions in figure EJa) could be interpreted as equal- 
magnitude high-speed regions whose strengths have been smeared by other less-dominant 
large-scale dispersive motions in the averaging process. Thus, a detailed description of 
physical processes far away from ajx is difficult to ascertain. One may be tempted to 
believe from the conditional average fields that t hese structure s are a ligned in the stream- 



wise direction; instantaneous visualisations, see iMontv et al\ ()2007() . suggest that these 
are in fact meandering structures, which would have been obscured in the averaging 
over the periodic domain in our study. However this apparent meandering coherence 
could equally well be interpreted as adjoined regions of shorter coherence, which are in- 
dividually well-captured by the conditional averaging. The spanwise scale appears to be 
approximately d. 

Figure [11^6) shows the spanwise structure of the relationship between ul and us- near 
the wall, Ul and us are positively correlated, but above z/S w 0.1, they are negatively 
correlated. This crossing point is different from z/S ~ 0.2 seen in figureSlJa). The discrep- 
ancy is resolved by noting the inclusive {Ai = ul — U > 0) and non-collocated (two-point) 
nature of the conditioning used for the averages {ul\Ai) and (wsjAi), as well as the single 
plane rather than integral representation. In contrast, the correlation i?, see (|3.8p . is con- 
structed from the one-point collocated statistic, {ui/us) = J {ul\ul){us\ul)p{ul) du^, 
which is not the same as J2Ai^L\A){us\A)P{A}. The opposite u^-us configuration of 
the weaker flanking regions in figure [5j&) suggests that these too experience phase rever- 
sal. 

Figure [5][c, d) shows the streamwise-wall- normal structure of the dominant large-scale 
motion and the relationship between the small and large scales. Despite a filter size of d, 
the coherence indicated in the figure suggests a wavelength of order 6(5, suggesting that the 
very long structures are the dominant contributors to the ul{xx) ~U < signal. Clearly 
the large-scale coherence has a streamwise phase that is dependent on the wall-normal 
location, at least where the coherence is strongest. Close to the wall this phase variation 
is weak, while the conditional averages with a;^ = (0, 0, 1000) 4^ Xx/6 ~ (0,0, 0.5) shown 
in figure[6]show that far from the wall the phase variation with increasing z/5 is also weak, 
but close to 7T rather than zero. In the intermediate region, the phase changes rapidly 
with wall-normal distance. The contours of constant \ug — {ug)\ in figure [5jrf) reveal a 
surface that cuts through the large-scale, low-speed region at a diagonal such that the 
region of negative correlation is larger where Ax < but smaller where Ax > 0. The 
angle of this separatrix, at least in the af orementioned interin ediate region, is about 14°, 
suggesting that the modulation reversal ( Mathis et al. 20091 ) is related to the structure 



inclination angle CMarusic fc Heueij|2007l ) 



Figure [5jc, d) and PjAi} w 1/2 suggest the stylised picture of a streamwise train of 
alternating high-speed and low-speed regions with the shape in figure Eljc). We propose 
that this sign change of ul — f/ determines the shape of the us region over a range 
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Figure 6. Conditionally averaged streamwise velocity fields in the (a) spanwise-wall-nor- 
mal and (6) streamwise-wall-normal planes, with {u'^\A2){Ax) and {u'^\A2){Ax), where 



A2 is the low-speed event ul{xx) — U < 
Xx/S = (0,0,0.5): filled contours of \{u+\A2) 



with probability P{A2} = 0.47 at 
- U+\ = 0.1,0.3,0.5,0.7, line contours of 



|(uj|^2) — (^'s)! ~ 0.01,0.02,0.03,0.04. Darker to lighter shades and dashed lines for nega- 
tive values. Vectors represent in-plane velocity components. Data from snapshot of Rct = 2 k 
channel flow LES, run Gib (table [1]). 



of wall-normal distances. Con sider the governing; equati on for m|, vifhich contains the 
production term —2ugduL/dx ( Reynolds fc Hussainlll972 ). Now, duh/dx < in between 
a high-speed region placed upstream {—Ax) of a low-speed region. Then, the production 
—2u^duL/dx > increases Ug, resulting in the picture [SJc?). The opposite mechanism 
applies in between a low-speed region placed upstream of a high-speed region, in which 
case —2ugduL/dx < 0, a backscatter of small-scale streamwise energy. This results in 
the quarter-phase shift between the and us region. 

A clearer picture emerges when we compute averages conditioned on the > to 
Ul < boundary, signalled by the event A^, = {duL/dx){xx) < (figure [71), interpreted 
as a quarter-phase streamwise shift of figure [51Jc, d). As expected, us is lowest precisely 
where Oul/Ox is minimum (at x^), figure 6). The bulge-like shape of the ul regions 
is preserved, see figure [TJ a), although with smaller sizes. Like the average conditioned 
on Ai, this figure is also reversible, with PjAa} = 0.48. 

The interactions further from the wall can be investigated by repeating the conditional 
averaging at x^/S = (0, 0, 0.5). A different relationship between the large and small scale 
activity emerges, as shown in figure [51 Instead of the 0-7t/2 phase difference close to the 
wall, |u J — U'^ I and |uj — {ug) \ are substantially out of phase, that is the phase difference 
is approximately n. This variation is in good agreeme nt with boundary layer resu lt s over 
a Reynolds number r ange of three decad es, namely [Bandvopadhvav fc Hussain ( 1984 ) 
and the recent work of Guala et al. ( 2009I) . From these other works, it would be expected 



that the large and small scales have close to zero phase difference very close to the wall; 
the LES formulation prevents us from confirming this point. 

A visual inspection of figure [5l gives an estimate that is of the order (0.8, 0.2, 0.05). 
Thus, the shear stress carried directly by these large-scale structures is estimated from 
this conditionally averaged picture as —{ulWl)^ ~ 0.04, and the large-scale streamwise 
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Figure 7. Conditionally averaged streamwise velocity fields in the (a) spanwise-wall-normal 
and (6) streamwise-wall-normal planes, with {u^\A3){Ax) and {ug\A3){Ax), where is 
the streamwise high-speed-to-low-speed boundary event {duL/dx){xx) < with probability 
P{A3} = 0.48 at Xx/S = (0,0,0.090): filled contours of |(u+|^3) - U+\ = 0.1,0.3,0.5,0.7, 
line contours of |(uj|y43) — {u'^)\ — 0.01,0.02,0.03,0.04. Darker to lighter shades and dashed 
lines for negative values. Vectors represent in-plane velocity components. Data from snapshot 
of Rct = 2k channel flow LES, run Gib (table[l]). 
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Figure 8. Conditionally averaged streamwise velocity fields in the (a) spanwise-wall-normal 
and (b) streamwise-wall-normal planes, with {u'^\A4){Ax) and (uJ|^4)(Aa;), where Ai is 
the streamwise high-speed-to- low-speed boundary event {duL/dx){xx) < with probability 
P{A4} = 0.48 at Xx/S = (0,0,0.5): filled contours of \{ut\Ai) - U+\ = 0.1,0.3,0.5,0.7, line 
contours of |(mJ|j44) — (ug)! — 0.01,0.02,0.03,0.04. Darker to lighter shades and dashed lines 
for negative values. Vectors represent in-plane velocity components. Data from snapshot of 
Re-r = 2k channel flow LES, run Gib (table[T|. 
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Figure 9. Spanwise and wall-normal velocities, and Reynolds stress associated with the con- 
ditionally averaged streamwise velocity in the (a) spanwise-wall-normal and {b, c) stream- 
wise-wall-normal planes, with {u'^\Ai){Ax) and {u'^\Ai){Ax), where Ai is the low-speed event 
ULixx) ~U < with probability P{Ai} = 0.53 at x,,/6 = (0,0,0.090): (a) filled contours 
of \{vt\Ai)\ = 0.04,0.08,0.12,0.16; (6) filled contours of \{w+\Ai)\ = 0.03,0.06,0.09,0.12; (c) 
filled contours of \{w+\Ai){{u+\Ai) - U+)\ = 0.02,0.04,0.06,0.08; and (a-c) line contours of 
|(Mg^|j4i) — (Mg)| = 0.01. Darker to lighter shades and dashed lines for negative values. Vectors 
represent in-plane velocity components. Data from snapshot of Rct = 2 k channel flow LES, run 
Gib (table [3. 



intensity is estimated as ~ 0.64. The relative magnitude and phase of the other 

velocity components associated with the large scale structure can be confirmed by looking 
at equivalent conditional averages for the other velocity components. Note that it is 
extremely difficult to obtain this sort of data experimentally, so we are effectively using 
the LES data in a predictive capacity to complete the description of the trends in the 
structure of the large-scale motions. 

Figures [51-[T1T a. b) reveal that while the spanwise velocity exhibits a footprint consis- 
tent with the implied swirl of figure [5J such that vl = statistically on the conditioning 
plane, but with a wall-normal phase variation similar to that exhibited by the stream- 
wise velocity u^. By contrast, the wall-normal velocity velocities in figures IMTTTc) are 
substantially in phase in the wall-normal direction, independent of the conditioning cri- 
terion. 

The conditional averaging technique can be taken one step further to demonstrate that, 
despite a large wall-normal footprint in the streamwise velocity, the large scales will con- 
tribute locally to the mean shear stress because of the relative wall-normal phases of the 
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Figure 10. Spanwise and wall- normal velocities, and Reynolds stress associated with the con- 
ditionally averaged streamwise velocity in the (a) spanwise-wall-normal and {b, c) stream- 
wise-wall-normal planes, with {u'^\A2){Ax) and {u'^\A2){Ax), where A2 is the low-speed event 
ul{xx) - U < with probability ^{^2} = 0.47 at Xx/S = (0,0,0.5): (a) filled contours 
of \{vt\^i)\ = 0.04,0.08,0.12,0.16; (b) filled contours of \{w+\Ai)\ = 0.03,0.06,0.09,0.12; (c) 
filled contours of \{w+\Ai){{u+\Ai) - U+)\ = 0.02,0.04,0.06,0.08; and (a-c) line contours of 
|(Mg^|j4i) — (Mg)| = 0.01. Darker to lighter shades and dashed lines for negative values. Vectors 
represent in-plane velocity components. Data from snapshot of Rct = 2 k channel flow LES, run 
Gib (table [3. 



large-scale streamwise and wall-normal velocities. Figures I MTTT c) show the conditioned 
shear stress distributions associated with the large scales, that is the product of the con- 
ditional averages of and w^. A conditionally averaged uj^wl yields a different answer 
that includes the effects of smaller structures since {ulWl\Ai) ^ {ui\A]) i^wi\A]) . The 
relative phases of and wj^ imply that there will be a contribution to the mean shear 
stress close to the wall. Note that this result is subject to the success of the conditional 
averaging in capturing streamwise coherence in and wl, but the trend is believed to 
be robust. 



4. Summary and Conclusions 

We have designed a series of LES runs that are well-suited to the investigation of large- 
scale structures in a long channel. The observations from this study lend themselves to 
an interpretation of the ubiquitous influence of large-scale structures in turbulent chan- 
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Figure 11. Spanwise and wall-normal velocities, and Reynolds stress associated with 
the conditionally averaged streamwise velocity in the (a) spanwise-wall-normal and 
(b, c) streamwise-wall-normal planes, with {u'j^\A4){Ax) and (uJ|j44)(Aa;), where 
A4 is the streamwise high-speed-to-low-speed boundary event {duL/dx){xx) < 
with probability P{^4} = 0.48 at x^/S = (0,0,0.5): (a) filled contours of 



Kw+|^i>| = 0.04,0.08,0.12,0.16; (b) 
filled contours of \{wj^\Ai){{uf^\Ai) — 



filled contours of \{iS+\Ai}\ = 0.03,0.06,0.09,0.12; (c) 
U+)\ = 0.02,0.04,0.06,0.08; and (tt-c) line contours of 
l(ujl^i) — — 0.01. Darker to lighter shades and dashed lines for negative values. Vectors 

represent in-plane velocity components. Data from snapshot of Rcr — 2k channel flow LES, run 
Gib (table 0. 



nel flow, consistent with th e boundary layer ex periments of iBandvopadhvay &: Hussain 
(|l984[) : [Ma"this et all (|2009f ) : iGuala et all (|2009[) . 

Using simultaneous time and spatial data to construct the spatio-temporal spectrum, 
we compute the dispersion relation, and show that departure from Taylor's frozen- 
turbulence hypothesis is noticeable near the wall, z/6 < 0.08, and for large scales, > 6S, 
in channel flow at Rcr = 2 k. This is consistent with the footprint of the very large scales 
reaching down to the wall. The opposite effect — that the large scales convect slower than 
the local mean — is also observed for these large scales away from the wall {z/6 ~ 0.5). 
This, too, is consistent with the conditional averages that show that these large-scale 
structures reach far away from the wall. 

Flow fields constructed from conditional averages confirm the extent of the influence 
of scales with A^: w 66 and reveal the wall-normal dependence of the spatial relationship 
between the high-speed large-scale region and the underlying small-scale intensity. The 
LES results appear to underline that the apparent amplitude modulation effect is better 
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described in terms of the spatial phase between the large and the small scales; indeed our 
correlation coefficients, Rp and Rt, are formulated in terms of this phase relationship. In 
this context, the zero in the correlation between large and small scale activity described 
in ii l3.2l can be interpreted as the location where, on average, the corresponding signals are 
7t/2 out of phase. A picture emerges in which the small scales are nominally in phase with 
the large scales near the wall and tt out of phase further from the wall. In the intervening 
region, where the phase difference is approximately 7t/2, the small scales track the sign of 
dui^/dx, per the small-scale production/backscatter term —2u^duL/dx. We can interpret 
the small scale activity, us, in terms of local structure, noting that the inferred locus 
of the maximum small scale energy corresponds to the well-known structure inclination 
angle of approximately 10-20°, at least at this Reynolds number. 

We conclude by observing that the very large scales appear to dictate some of the 
turbulence behavior close to the wall, with interesting implications for the "top-down 
versus bottom-up" debate concerning the dynamical significance of different regions of 
the flow. It is perhaps more correct to say simply that the two regions can now be 
determined to be inextricably linked. 
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